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ABSTRACT 


For Terahertz (THz) imaging to be useful outside of a laboratory setting, inexpensive, yet 
sensitive, detectors such as uncooled microbolometers will be required. Metamaterials 
can improve THz absorption without significantly increasing the thermal mass or using 
exotic materials because their absorption is primarily dependent on the geometry of the 
materials and not their individual optical properties. Finite Element (FE) simulations 
revealed that an array of squares above a ground plane separated by a dielectric is 
efficient, yet thin. Metamaterials were fabricated and their absorption characteristics were 
measured using a Eourier Transform Infrared Spectrometer indicating that the EE 
simulations are accurate. Metamaterial structures tuned to a quantum cascade laser (QCE) 
illuminator were incorporated into a bi-material sensor, which was used for detection of 
THz radiation from the QCE source with good sensitivity. In the case of 
microbolometers, a bolometric layer needs to be embedded in the metamaterial to form a 
thin microbridge. Simulations indicated that if the bolometric layer was resistive enough 
or close enough to the ground plane, then absorption would be largely unaltered. 
Metamaterials with a conductive Titanium layer embedded into the dielectric spacer were 
fabricated and measured with an ETIR, confirming this behavior. 
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I. 


INTRODUCTION 


A. TERAHERTZ IMAGING 

In recent years, there has been a surge of renewed interest in the terahertz (THz), 
or Far-Infrared (Far-IR), spectral range. The THz band (Figure 1), so named because the 
frequencies involved are typically between 0.1-10 x lO'^ Hz (3 mm-30 pm wavelength), 
has some interesting real world imaging applications. Terahertz easily penetrates dry, 
non-polar, low-conductivity materials (such as common packaging materials and 
clothing), but is readily reflected by metallic objects and strongly absorbed by water. This 
limits the effective range, especially above 1 THz, to short-range applications. 
Additionally, THz radiation is non-ionizing and does not penetrate deeply into the skin, 
making it potentially safer than X-ray imaging technologies and much less likely to 
damage objects [1]. 
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Figure 1. Electromagnetic spectrum showing the location of the terahertz band. 

From [2] 


Therefore, two of the main interests for THz imaging is security screening [3] and 
medical imaging [4], [5]. In the case of security screening, metallic objects such as knives 
are easily observed and explosives are known to have unique absorption spectra in this 
region allowing them to be identified inside packages without damaging the contents. In 
the case of medical imaging, much of the interest lies in the absorption of THz by water. 
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Skin cancers tend to have more blood vessels than surrounding tissue and thus more 
water content. This can allow one to image eancers with THz [4], [5], perhaps even in 
real-time during surgery. Cultural heritage conservation [6] is another fascinating 
application, where THz imaging ean be used to identify art obscured by other materials 
(sueh as multiple layers of paintings) or help validate the authenticity of priceless objeets 
without damaging them. Additionally, quality eontrol of food [7] can make use of the 
same prineiples as medical imaging. Food can be easily imaged inside packages to track 
how well it has dried, or even in real time to traek the effieieney of the proeess [7]. 
Unfortunately, detection in the THz range is quite ehallenging due to the low photon 
energy necessitating excessive cooling for optical schemes and the high frequeney 
signifieantly reducing the effeetiveness of electronic schemes; these issues ereate the 
well-known “Terahertz gap” [8]. Sensors are available in this area, but are prohibitively 
expensive and bulky, largely restricting them to the researeh or remote sensing domains. 

B. TERAHERTZ IMAGING USING UNCOOLED THERMAL DETECTORS 

Uncooled thermal detectors are a well-established, eompact and inexpensive 
technology for imaging in the infrared, especially in the form of mierobolometer focal 
plane arrays (FPAs). Indeed, these sensors naturally evolved in snakes in the form of pit 
organs leading to a high sensitivity, albeit low resolution, infrared imaging system for 8- 
12 pm radiation and to a lesser extent the 3-5 pm range [9]. As this sensor is foreed to 
operate at the temperature of the snake (whieh is the same as the thermal baekground 
since they are cold-blooded), it is not surprising that a similar uncooled sensor would 
perform well when engineered by humans faeing similar limitations. 

All thermal sensors, ineluding those in snakes, work through a thermally isolated 
sensitive area heating up when their thermodynamic equilibrium is disrupted by a change 
in ineident eleetromagnetic radiation. Thermal sensors differ in the way they measure this 
temperature change; mierobolometers, one of the more dominant imaging technologies 
for uncooled 8-12 pm imaging, work using a layer of material in the thermally isolated 
region (microbridge) that experienees changes in electrical resistance with temperature. 
While all materials ehange resistanee in response to a temperature shift, for eertain 
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materials such as titanium (Ti) and Vanadium Oxides (VaOx), the magnitude of this 
resistance change is much greater making them especially good choices [10]. 

As demonstrated in Figure 2, there is a large amount of power in the 8-12 pm 
infrared range at common temperatures, but power in the TFlz region is comparatively 
low. Therefore, active sources are commonly used as illuminators [3]-[7], [11]-[13] for 
THz imaging systems. Possible active sources include free electron lasers (FELs) [14], 
quantum cascade lasers (QCLs) [15]-[17], diode based oscillators [18], and nano¬ 
electrodes in a photoconductive photomixer [19]. QCLs, in particular, are both compact 
and powerful, providing a useful narrowband source for imaging. One drawback is that 
QCLs operate at cryogenic temperatures. 

Wavelength (nm) 



Figure 2. Blackbody distributions for various temperatures with the TFlz band shaded in 
blue and IR in red. THz makes up only a small fraction of blackbody 
emission at common temperatures. From [20]. 

THz images have been successfully produced using a microbolometer camera and 
QCL source in a laboratory setting [11]-[13], allowing one to image concealed objects, 
as shown in Figure 3(a). However, the sensitivity of these microbolometer cameras is 
limited because the silicon nitride (SiNx) microbridge is not optimized for high THz 
absorption [21], [22] as illustrated by the absorption curve in Figure 3(b). Indeed, 
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uncooled IR microbolometers generally have a noise equivalent power (NEP) in the THz 
range of approximately 300 pW/Hz'^" [21] compared to an NEP of 14 pW/Hz'^" for the 

IR range [10] implying low responsivity at THz frequencies. Therefore, improving the 
absorptivity of these microbolometer pixels should improve their performance 
considerably. Metamaterials, with their subwavelength features, can possess absorption 
characteristics not found in their constituent materials, which allows fabrication of 
structures with optimized absorption at THz frequencies. These structures can be 
incorporated in thermal sensors for enhancing the sensitivity at THz frequencies. 




Figure 3. Knife covered in opaque white plastic tape (a) is visible in THz. 

Absorption (b) of a 1 pm thick Silicon Nitride membrane, showing poor 
absorption in THz. (b) From [20]. 


C. OVERVIEW OF PRESENT RESEARCH 

In the following chapters, I present the research progression from designing, 
fabricating, and measuring metamaterial absorbers alone, integrated into bi-material 
sensors, and with a bolometric layer added [23]-[26]. Although each chapter is largely 
based off a single published work, some of the material has been reorganized for a more 
logical presentation. Each phase of the research required Finite Element (EE) modeling, 
clean room fabrication (performed by my colleague Dr. Fabio Alves), and finally 
measurements in the laboratory using a Quantum Cascade Laser (QCL) source and a 
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Fourier Transform Infrared Speetrometer (FTIR). The results of this researeh demonstrate 
that metamaterials are a promising ehoiee for improving the THz absorption of thermal 
sensors without eompromising their real-time imaging performanee. They also offer 
guidanee on how to properly integrate them into mierobolometer pixels, allowing one to 
eonstruet highly sensitive THz mierobolometer foeal plane arrays in the future. 

In Chapter II, metamaterials eonstrueted from fabrieation-friendly materials are 
investigated and optimized. FE analysis is used to guide the fabrieation of metamaterial 
absorbers optimized to absorb nearly 100% of ineoming THz power at a speeifie resonant 
frequeney. The metamaterials are eonstrueted from A1 and Si-rieh SiOx; both of these 
materials are very suitable for eonstruetion of sensors, espeeially bi-material sensors as 
elaborated on in Chapter III. The speeular refleetion of THz radiation by these 
metamaterials is measured using our FTIR and used to determine the THz absorption of 
these layers (transmission and higher order seattering are negligible). These results 
provide us with the information needed to tune metamaterial absorbing layers to the 
desired frequeney with nearly 100% absorption. Finally, the emission properties are also 
investigated showing that these materials obey Kirehoffs law. The majority of the 
eontent of this ehapter was published in largely the same form in Optical Engineering 
[23], and presented at the international soeiety for op ties and photonies (SPIE) Defense 
and Seeurity Sereening eonferenee in April 2012. The emission seetion was published in 
essentially the same form in Optics Express [24], and presented at the SPIE Photonies 
West eonferenee in Eebruary 2013. 

Chapter III uses the knowledge gained in Chapter II to design and integrate 
metamaterial absorbers into bi-material pixels, hirst, the basie physieal limitations of 
thermal sensors are explored and the need for layers with a small thermal eapaeitanee is 
established. This seetion is extended from a paper published in Optical Materials Express 
[25]. After this, the details behind design, fabrieation, and measurement of bi-material 
pixels integrated with THz metamaterial absorbers is explored. EE modeling is used to 
optimize the thermomeehanieal sensitivity of the pixels, while the data from Chapter II is 
used to optimize absorption. The sensors are fabrieated using standard eleanroom 
mierofabrieation teehniques and their performanee is eharaeterized using a QCE souree 
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combined with an optical readout setup to determine the deflection under THz 
absorption. This section demonstrates that the metamaterials in Chapter II can be 
integrated into bi-material pixels, increasing their THz absorptivity. Other than the first 
section, this work is contained entirely within a paper published in Optics Express [26], 
and was also presented at the SPIE Photonics West conference in February 2013. 

Although bi-material pixels face the same thermal challenges that microbolometer 
pixels do, suggesting that metamaterials will enhance their performance, they do not have 
the same structure. Indeed, the effect of the bolometric layer in a microbolometer pixel 
must be explored to ensure that metamaterials can incorporate this layer without having 
their absorption adversely affected, or the thermal properties of the pixels compromised. 
Therefore, to investigate this metamaterials with an embedded bolometric layer (in this 
case, Ti) were fabricated. These metamaterials were fabricated with the bolometric layer 
at several different depths and the subsequent FTIR measurements demonstrated that if 
the bolometric layer is close to the ground plane, the absorption is largely unaffected. 
This information was used to validate the FF model, which also suggests that a 
bolometric layer with a low sheet resistance will also have a smaller impact on the 
maximum absorption amplitude. Finally, different possible configurations are explored 
for their various advantages and disadvantages and future work in the area is suggested. 
Other than the last sections, this chapter is published practically unaltered in Optical 
Materials Express [25]. 

The progression from metamaterial absorber design, integration into thermal 
sensors, and analysis of the impact of a bolometric layer represent significant progress 
towards the successful integration of THz metamaterial absorbers to enhance the 
performance of thermal sensors, specifically microbolometers. Chapter II provides the 
general information necessary when constructing these layers. The results from Chapter 
III confirm that these metamaterials can be integrated into thermal sensors and that 
constituent materials that maximize the thermal sensitivity of the sensor can be used 
while still maintaining high THz absorptivity and a small thickness, highlighting the 
primary advantage of metamaterials. Finally, Chapter IV confirms that a bolometric layer 
can be integrated and suggests that this layer should be placed as close as possible to the 
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ground plane and have as small of a sheet resistanee as possible. Ideally, these advanees 
will lead to highly sensitive THz mierobolometer eameras using metamaterial absorbers 
for praetieal real-time imaging applieations. 
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II. METAMATERIAL ABSORBERS 


Note: The majority of the content of this chapter was published in largely the 
same form in Optical Engineering [23], and presented at the international society for 
optics and photonics (SPIE) Defense and Security Screening Conference in April 2012. 
The emission section was published in essentially the same form in Optics Express [24], 
and presented at the SPIE Photonics West Conference in February 2013. 

Metamaterials are artificial structures with subwavelength features that can 
exhibit overall optical properties not characteristic of their constituent materials. Indeed, 
they can even have characteristics that are not found in nature, such as a negative 
refractive index [27]. Absorption in these structures, especially through resistive losses, 
can be a vexing problem for applications such as superlenses [28]. However, certain 
applications can benefit greatly from enhanced absorption, such as thermal sensors. 
Eandy et al. [29] initially envisioned a material which would both not transmit light yet 
would be impedance matched to free space through adjustment of the dielectric and 
magnetic response of the metamaterial. This would result in a material that could not 
reflect or transmit light, and therefore must absorb all incident radiation. In practice, this 
typically is only achievable over certain frequency bands, depending on the geometry of 
the subwavelength features. Materials with this property have been realized across the 
electromagnetic spectrum [30]-[43]. In principle, one can fabricate a metamaterial for 
any frequency within this range of interest by an appropriate adjustment to the geometry 
of the subwavelength features. Eortunately, this tuning is only weakly dependent on the 
materials used, so these metamaterial absorbing films can be constructed using standard, 
mature MEMS microfabrication techniques. These metamaterial films typically consist of 
a periodic-patterned metal layer, a dielectric spacer, and a ground plane [30]-[43]. In this 
work, a periodic array of metal squares was used as shown in Eigure 4. Metal squares 
were chosen due to the ease of their fabrication and thin dielectric spacer requirements. 
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A. 


SIMULATION 


While empirical studies are, of course, essential to evaluate the properties of 
metamaterials, numerical modeling can very effectively predict the properties of a 
metamaterial absorber without the need to fabricate every possible geometry. Reflection 
is suppressed as incident EM waves excite localized surface plasmon resonances, leading 
to resistive losses in the metal and additional losses in the dielectric layers converting 
some of the incident energy to heat, allowing the incident radiation to be absorbed [30], 
[32], Resonant absorption frequency (v) and absorption magnitude can be understood 
using explanations such as impedance-matching to free space [29], an equivalent resistor 
inductor capacitor (RLC) circuit [30], [32], confined TM cavity modes [44], interference 
of multiple reflections [45], and transmission lines [46]. However, due to the complex 
nature of metamaterials, these models are typically guided by finite element (FE) 
modeling (of at least the patterned metallic layer) [30], [32], [44]-[46]. Therefore, 
understanding numerical models of metamaterial behavior is essential to predicting their 
properties. In this work, I used COMSOE Finite Element modeling software. COMSOE 
is a multiphysics package with many different modules available, but the one used for 
analyzing metamaterials is the radio frequency (RE) module. Great care must be taken 
when constructing 3D models, so that one can preserve the physics while still obtaining a 
result in a reasonable amount of time. 

When creating a FE model, boundary conditions and meshing are crucial 
considerations. The geometry itself is usually fairly straightforward to create and so is of 
little concern. I should emphasize that the geometry should match the fabricated structure 
as accurately as possible, as even minor changes in geometry can lead to variations in 
absorption characteristics. A typical rule of thumb for FE analysis is that elements should 
be no further than k/5 apart, which can result in an excessive amount of mesh if a very 
large volume of metamaterial is simulated in a 3D model. To simplify the situation, we 
can exploit the fact that metamaterials are periodic structures and thus any 
electromagnetic solution must satisfy periodic boundary conditions. Therefore, only a 
single unit cell of metamaterial (see Figure 4(a)), with appropriate periodic boundary 
conditions, is necessary for FE analysis. For radiation normally incident to the surface, 
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the electric and magnetic fields of an incident plane wave can be broken into two 
polarizations with fields normal to the boundaries. In this case, a perfect electric 
conductor (PEC) will enforce periodic boundary conditions for electric fields normal to 
the boundary, while a perfect magnetic conductor (PMC) will do the same for magnetic 
fields normal to a boundary. In the case of radiation not normally incident to a surface, 
Floquet boundary conditions are appropriate as they attach the necessary phase factor to 
the radiation [47]. However, in most thermal sensor applications radiation is near normal 
incidence making normally incident radiation a good approximation. 




Scattering Boundary Condition 

(a) (b) (c) 

Figure 4. Unit cell (a) of a metamaterial with a single absorption peak. The model also 
incorporates the unintentional over-etch into the SiOx during the fabrication. 
Schematic representation (b) of the unit cell in the FE modeling program with 
boundary conditions and THz plane wave at normal incidence. Time- 
averaged power flow (c) near the resonant absorption frequency and arrow 
lengths relate to magnitude in log scale. From [23] 


The model must also allow electromagnetic radiation to interact with the 
metamaterial absorbers. To accomplish this, one must choose the proper boundary 
conditions for the surfaces intended to allow radiation to penetrate through. One surface 
must somehow generate, or allow penetration of, incident radiation and also allow 
reflected radiation to leave the model or absorb it without allowing any radiation to 
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reflect back in. The second surface, typically “below” the ground plane, must allow any 
transmitted radiation to leave the model or absorb it, also without any radiation reflecting 
back into the model. Several popular methods available in COMSOL are scattering 
boundary conditions, internal ports, and perfectly matched layers (PMLs). Scattering 
boundary conditions are meant to approximate an infinite space or “open” boundary for 
radiation to propagate out into. They work as first order absorbing boundary conditions 
for plane, cylindrical, or spherical waves. Port boundary conditions are similar in that 
they are perfectly absorbing boundary conditions for general modes of known shape. 
Both of these boundaries can also allow the same types of waves to penetrate the 
boundary into the model, supplying the incident radiation. The power is easily controlled 

by adjusting the Poynting vector and can be calculated using A[S-n^ for plane 

boundaries. While not exactly a boundary condition, a PML can also be used to 
approximate an open boundary. A PML is actually very similar to the initial idea of a 
metamaterial absorber; it is impedance matched to the material that it is adjacent to so 
that it can eliminate reflection and uses a complex coordinate transformation to cause an 
incident electromagnetic wave to decay as it propagates resulting in essentially total 
absorption [47]. The PML actually has a finite thickness and must have additional mesh 
elements, and should usually be coupled with a scattering boundary condition to 
minimize back reflections. All of these tools tend to work very well for a particular wave 
mode and propagation vector and then start to allow some reflection back into the model, 
so care should be taken in their implementation. A schematic of this setup is shown in 
Figure 4(b). 


Once boundary conditions, geometry, and meshing are finished, the wave 
equations for the model to solve must be specified. The default option for COMSOL is 
“relative permittivity” which solves the EM wave equation of the form; 
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Thus, this equation admits relative permeability and permittivity as well as 


conductivity. This is convenient for materials that can be defined entirely by permeability 
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and permittivity, such as free space, or those where conductivity dominates, such as most 
metals in THz. For some metals, frequency-dependent conductivity is needed if it varies 
appreciably within the frequency range of interest (such as Au [44]). In the case of A1 
used in our metamaterial structures, the relaxation is fast enough that the Drude model 
approximates to the DC conductivity [49]. It should be noted that thin films of metals do 
not have the same conductivity as bulk metals and also may have variable conductivities 
depending on deposition technique and surface [49]. COMSOL does have a “Drude- 
Lorentz dispersion model” option, but this approximation renders it unnecessary for A1 as 
absorption remains largely unchanged. The second convenient setting is “Refractive 
index” which solves the EM wave equation of the form; 

V xE^-kls^E = Q (2) 

Clearly, this form assumes relative permeability is 1 and conductivity is contained 
in the relative permittivity via the extinction coefficient [46]. It then calculates the 

permittivity using the complex refractive index: s^={n- ik )^, using the usual 

convention that a positive “k” refers to loss while a negative “k” is gain. This form is 
useful for dielectrics, where the permeability can largely be ignored, and the literature 
commonly specifies the complex refractive index. 

Once the model construction is complete, the solvers must be chosen. The models 
in this work use a frequency domain solver, which simply solves the EM wave equations 
at specified frequencies. This results in the E and B-field solutions allowing one to 
determine information such as power flow (see Eigure 4(c)). More resolution simply 
correlates to a longer computation time, but more precise absorption curves. Of greater 
consequence is the type of solver used. Direct solvers are preferable if the modeling 
computer possesses sufficient memory as there are no convergence issues, yet indirect 
solvers can solve denser meshes with the same resources but can have convergence 
issues. Thus, the choice of solver and settings will largely be determined on a case-by- 
case basis. Direct solvers were used for the models in this work, due to the computational 
resources available. In addition, parametric sweeps can be used to run several different 
configurations of metamaterials together. 
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Finally, once the model has finished its run, one must extraet absorption 
eharaeteristies. In all oases, when magnetio losses are largely ignorable in the oonstituent 
materials, a oonvenient way to extract absorption is through performing volume 
integrations of the resistive losses in the metamaterial. Non-linear effeots are largely 
ignored in these models, so the inooming power oan be set to 1 W and the absorption can 
simply be read off as the Watts of power absorbed in the metamaterial. This also has the 
advantage of identifying which layers are aotually heating up, allowing one to more 
effeetively model metamaterials that are not thermally linked as in Figure 5. 



Figure 5. Modeled absorption eharaeteristies of the three layers of a metamaterial with 

single absorption peak. From [23] 


Additionally, when using ports, one can also extract transmission and reflection 
coefficients by inverting the S parameters. This method is also important for the 
refleetion model that relies on the resonator layer having some frequency dependent 
refleetion eoeffieient [45]. This, of eourse, assumes that higher order seattering is largely 
ignorable whieh should be true for sub-wavelength features of this kind and is known to 
be small for metamaterials in other frequeney ranges [46]. In praetiee, transmission 
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through a metamaterial is vanishingly small due to the optieally thiek ground plane of 
approximately 100 nm. The skin depth of A1 layers is approximately 80 nm for 1 THz, 
and both models and measurements show no signifieant transmission. Finally, refleetion 
and transmission eoeffieient teehniques and resistive loss ean be used together to validate 
that models are working eorreetly [30]. 

B. FABRICATION 

Metamaterial absorbers, whether split ring resonators or arrays of squares, 
typieally eonsist of a patterned metal layer, a dieleetrie spaeer of some kind, and a ground 
plane. Variations ineluding a bolometrie layer and integration into bi-material sensors 
will be diseussed in later ehapters. Fabrieation of these layers is easily done using 
standard mierofabrieation teehniques. This was performed at Oak Ridge National 
Laboratory by Doetor Fabio Alves and Assistant Professor Dragoslav Grbovie. At every 
point in the fabrieation, the materials were earefully monitored for relevant parameters 
sueh as residual stress and DC eonduetivity. These measurements were performed with 
stylus profilometry and optieal interferometry in the ease of thiekness, four-point probe 
measurements in the ease of eonduetivity, and wafer eurvature measurements using a 
stylus or optieal profilometer to evaluate residual stress. As transmission is essentially 
zero due to the thiek ground plane of the metamaterial [30], [48], one does not need to 
use low eonduetivity substrates as with thin metal films [49]. If eonduetivity of the wafer 
is a eoneern, a layer of insulating dieleetrie ean be deposited. Therefore, a single side 
polished (SSP) Si wafer of any thiekness or doping eoneentration should be suitable, but 
thieker wafers are less fragile. Generally 500 pm thiek wafers were used due to the low 
priee and durability of these wafers. 

The first step is to deposit the metallie ground plane. We typieally used A1 
deposited using an e-beam evaporator due to the fast deposition rate, good adhesion, and 
good eonduetivity. Typieally, these A1 layers were measured to have a DC eonduetivity 
of 10^ Siemens/meter for layers of approximately 100 nm. The bulk eonduetivity of A1 is 
over 3 times larger than that of the A1 films, therefore measurement with a four-point 
probe is essential. The seeond step is the dieleetrie spaeer, whieh in this work was SiOx. 
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Specifically, this was a low stress oxide deposited with a Plasma Enhanced Chemical 
Vapor Deposition (PECVD) system at 350°C with more Si than the standard Si02 
formulation (-13 MPa compared to -140 MPa). Dielectric film parameters can vary 
between machines and even through small changes in machine parameters; therefore, 
great care was taken to reproduce the same SiOx for every phase of the project. Einally, a 
second A1 layer is deposited on top of this dielectric. The layer is deposited using exactly 
the same settings as the ground plane metal. To pattern this layer and form a 
metamaterial, it must be etched. Eirst, after the wafer is prepared, photoresist is spun on 
using a spin coater. After some baking steps specific to the photoresist, the wafer is 
placed into a contact aligner. A mask with the metamaterial absorber resonators is then 
placed on top of the wafer and the photoresist is exposed to UV light. After this, the 
photoresist is developed and processed further, again with steps specific to the type of 
photoresist. In this work, positive photoresist was used, so unexposed areas remained on 
the wafer while exposed areas were dissolved. To remove the unwanted, now exposed Al, 
an Ar plasma sputter etcher was used. The sputter etcher accelerates Ar ions at the wafer, 
removing both photoresist and exposed Al. This etcher is anisotropic (directional) which 
reduces the “undercut” one would get from isotropic etching techniques such as acid 
etching. The remaining photoresist protects the Al that we wish to stay on the wafer. 
Unfortunately, this process is not very selective between Al and the underlying dielectric 
and the etching must occur long enough to go all the way through the Al layer. This 
results in an effect known as “over-etch” where some of the dielectric is also removed as 
illustrated in Eigure 4(a). The effect of over-etch can be accounted for in simulations by 
simply adding to the geometry. Einally, the photoresist is removed using appropriate 
chemicals and an isotropic oxygen plasma “barrel” etcher. The result is a wafer with 25 
different metamaterials pictured in Eigure 6. 

It should be noted that other materials and processes can be used. Indeed, acid 
etching (aka “wet” etching) is much more selective and could result in no over-etch into 
the dielectric. If the etching is carefully controlled and reproducible, it may well be a 
better choice than sputter etching. Other suitable choices for the Al etching include 
reactive ion etching and liftoff. The advantage of reactive ion etching is its better 
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selectivity towards Al. In the case of liftoff, the photoresist is deposited first. Inverted 
masks are used such that the unwanted Al is deposited on the photoresist and is removed 
with it. This requires no metal etchant at all, and is therefore faster and easier, albeit with 
less precision. Also, other materials can be used for metals include Au or Ag or 
dielectrics such as SiNx and polyimides. 



Figure 6. One of the fabricated wafers with 25 different metamaterials arranged radially 
for automated reflection measurements using a PIKE MappIR FTIR 

accessory. From [23] 


C. MEASUREMENT 

To evaluate the THz optical properties of the metamaterial structures, a Fourier 
Transform Infrared Spectrometer (FTIR) was used. This device consists of a Michelson 
interferometer with a moving mirror in one arm; this results in periodic oscillations in 
amplitude whose frequency depends on the wavelength of radiation and velocity of the 
mirror. Thus, the FTIR’s thermal detector can measure every frequency simultaneously 
and extract the frequency data with a Discrete Fourier Transform (DFT). This powerful 
technique is applicable across the frequency spectrum, but is particularly effective for the 
Infrared region of the spectrum. The FTIR available at NPS is a Thermo-Nicolet Nexus 
870 fitted with a globar source (hot SiC filament) and a Deuterated Triglycine Sulfate 
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(DIGS) pyroelectric detector with a polyethelene window. This is the recommended 
configuration for Far-IR spectrometry (defined by Nicolet as 50-600 cm"' or 1.5-30 
THz); however, the 1.5-3 THz range is difficult to use due to the fact that the source 
power is relatively low as shown in Figure 7. For measurements in the Mid Infrared 
(Mid-IR) range (defined by Nicolet as 7400-350 cm'^ or 1.35-28.6 pm), the beam splitter 
is changed to KBr and the detector is changed to a Thermoelectrically cooled DIGS 
detector with a KBr window. For metamaterial absorption, the Far-IR configuration is of 
important as it contains the THz range of interest. In terms of FTIR settings and 
configurations, one must achieve a large enough signal to align the FTIR, yet a small 
enough signal to avoid saturation of the detector. In the default configuration of the FTIR, 
it can only measure transmission through a sample. 



Figure 7. Raw FTIR data for a metamaterial absorber, showing the excessive noise in 

the spectral region below 3 THz. 


Clever use of optics can allow one to measure reflection as shown in Figure 8; 
specifically, an FTIR accessory known as MappIR can be programmed to automatically 
measure either transmission or reflection at predefined points on a wafer resulting in a 


18 




significant time savings when dealing with a large number of metamaterials. In the case 
of emission measurements, the FTIR is programmed to see the metamaterials and a black 
carbon blackbody (prepared by depositing candle soot on a suitable substrate) as an 
external source. An off axis parabolic mirror collimates the thermal radiation from a 
small area of the sample on a hot plate and directs it into the FTIR. In all of these 
configurations, the THz must pass through the outside air, leading to frequeney 
dependent losses that cannot be ignored. The inner area of the FTIR is purged with elean, 
dry, low CO 2 air to preserve the optics and boost signal strength. Despite this a 
background spectrum must always be taken in FTIR measurements to subtract the non- 
uniform spectral characteristics of the system. Appropriate backgrounds depend on the 
measurement and configuration, but typically are air for transmission, thick uniform Au 
or A1 films for reflection, or black carbon blackbodies for emission. The FTIR can 
automatically compare these measurements to the sample or they can be done in the 
graphing software. Typically the output is given in “percent transmission” whieh 
essentially divides the sample points by the background points. This results in sizable 
errors when the signal is mueh smaller than the noise as seen for the 1.5-3 THz range in 
Figure 7. 



(a) (b) 


Figure 8. Schematic representation of the experimental configuration used for 

measuring the reflection (a) and transmission (b), before using the MappIR 
accessory. The accessory performs essentially the same measurements, but 
operates at an incident angle of 15 degrees. From [49]. 
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D. 


OPTIMIZATION AND TUNING OF SINGLE-BAND METAMATERIAL 
ABSORBERS 


The simplest metamaterial geometry contains only one type of square arranged in 
a periodic array. I refer to them as “single-band” absorbers because they have only one 
absorption peak in the frequency range of interest. As Figure 9 demonstrates, we see that 
the lowest order confined mode consists of a single oscillating dipole formed by the 
incident radiation interacting with the square with opposing currents and charges in the 
ground plane. 



Figure 9. FE model solution for 16 pm squares at the resonant absorption frequency. 

Electric (red) and magnetic (blue) fields for a metamaterial absorber at its 
resonant frequency. Note the opposing electric fields on either side of the 
square and the reversal of the magnetic field. 


This suggests that we should expect the wavelength of any trapped first-order 
modes should be proportional to the size of the structure in the direction of the incident 
E-field component parallel to the surface. Of course, this does not take into account the 
boundary conditions or the coupling between squares. Eigure 10(a) shows experimental 
results for single band absorbers with the same 21 pm periodicity, but incrementally 
increasing square size. As expected, the resonant frequency decreases as square size 
increases. Eigure 10(b) shows that the resonant frequency is roughly proportional to the 
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inverse of the square size, as expected due the fact that wavelength is inversely 
proportional to frequency. Therefore, modifying square size is an effective method for 
controlling the resonant absorption frequency. 




Frequency (THz) Inverse Square Size (nm ’) 



Figure 10. FTIR measurements (a) for varying square size, 21 pm pitch, and 1.47 pm 
SiOx . Linear fit (b) of inverse square size (1/s) to resonant absorption 
frequency (v). The effects of varying dielectric thickness (c) for 12 pm 
squares and 21 pm pitch and varying pitch (d) with 16 pm squares and 1.47 

pm SiOx. From [23] 


Now that the most effective frequency tuning method has been established, we must 
optimize the remaining geometric parameters to achieve the highest possible maximum 
absorption amplitude. Dielectric thickness and periodicity are known to be significant 
parameters for maximum absorption. As Figure 10(c) demonstrates, a higher dielectric 
thicknesses typically result in higher absorption. As no reduction in absorption amplitude 
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was observed for inereased thiekness in this range, it seems that all absorbers have an 
optimal thickness at or above 1.66 pm. Once past this optimal thickness, absorption 
magnitude will decrease with increasing thickness until another optimal thickness is 
approached. The reflection model of absorption suggests that, similar to thin film 
interference, there should be several optimal thicknesses. In the case of thermal sensor 
applications, we will see that in most cases these materials should be as thin as possible. 
Periodicity, or pitch, shows a similar relationship to dielectric thickness, showing more 
significant frequency shifts as demonstrated in Figure 10(d). Intuitively, a smaller pitch 
should correspond to a higher absorption due to a higher density of squares [30], [32] or 
better coupling to external radiation [44]. It should be noted that there is a relationship 
between pitch and dielectric thickness such that a given dielectric thickness will have an 
optimal pitch. As Figure 11 shows, simulations suggest that at certain thicknesses 
increasing pitch can increase absorption magnitude. While this is certainly relevant if one 
wishes to construct very thick absorbers, in the thin regime a decrease in pitch will 
always be desirable. Of course, any resonant absorption frequency shifts due to 
periodicity or dielectric thickness changes must be compensated for through tuning the 
square size. 

While metamaterial absorbers can be constructed using many different materials, 
the materials used do have a significant impact on the optimal dimensions for the 
metamaterial geometry. Metals, or equivalent high conductivity materials, must have 
sufficient thickness and conductivity to function as resonators. Additionally, some metals 
such as Au have conductivities that vary greatly across 1- 10 THz [44], [45], [48]. 
However, the A1 films used in this experiment have essentially the same conductivity 
across the 3-8 THz range [48]. The refractive index of the dielectric spacer is also 
significant, as it will relate to the capacitance of the squares or the wavelength of the 
trapped dipole modes, resulting in higher absorption frequencies for lower refractive 
indexes. Additionally, processing steps that over-etch into the dielectric cannot be 
ignored. This removes dielectric from around the squares, causing a shift in resonant 
absorption frequency. This could, of course, be used intentionally to shift resonant 
frequency and reduce the thermal mass of the absorbers if desired. 


22 




Figure 11. Simulations for 16 [xm square size and 21 qm pitch with two dielectric 

thicknesses showing that for thicker cases, increasing periodicity can, in fact, 

increase percent absorption. 


E. OTHER METAMATERIAL CONFIGURATIONS AND EMISSION 

CHARACTERISTICS 

While single band metamaterial absorbers are ideal for coupling with a single 
narrowband source such as a QCL for active THZ imaging, in some cases a broader 
absorption peak may be desired. Additionally one might require two absorption peaks, for 
instance when working with two separate QCL sources. Metamaterials can be fabricated 
with these properties by using more complex unit cells. A possible configuration is to 
simply alternate the squares resulting in a pattern reminiscent of a checkerboard, as 
illustrated in Figure 12(a) with squares dimensions of 13 and 17 pm and a dielectric layer 
thickness of 1.47 pm. Figure 12 shows the measured (b) and simulated (c) absorption of 
structure in (a) along with absorption spectra by removing one type of squares and 
keeping the rest of the parameters the same. 
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(a) 



Frequency (THz) 



(b) 



(c) 


« 1 ) 


Figure 12. Unit cell (a) of a dual band metamaterial configuration. Measured (b) and 
modeled (c) absorption for a dual-band absorber with squares that have a 4 
pm size difference compared to the same configuration with only one type of 
square. Comparison (d) of the absorption spectrum for dual band absorber 
with 2 pm square size difference, overall pitch of 42 pm and SiOx thickness 
1.47 pm to a single band absorber with an intermediate square size. From 

[23]. 


It can be easily seen from Figure 12(b) that each type of square is responsible for 
one of the two peaks of the combined structure absorption. For the configuration shown 
in Figure 12(a), the absorption measurement shows two well-separated absorption peaks 
at 4.25 and 5.5 THz similar to that obtained in the simulation. In addition. Figure 12(b) 
shows measured single-band absorption from two complementary structures fabricated 
with only one type of square while maintaining the same dimension and pitch. In this 
case, single-absorption peaks clearly match with the absorption peaks of the dual-band 
absorber, indicating that each set of squares in Figure 12(a) is responsible for generating 
the two peaks as predicted by the FE model (Figure 12(c)). As the difference in 
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dimensions of square size is reduced (16 and 18 qm), the absorption peaks begin to 
overlap and merge together as illustrated by measured data in Figure 12(d). Thus, by 
bringing the size of the squares closer together, it is possible to produce a broader 
absorption peak than a single square configuration of intermediate size (17 pm). In both 
the measurements and the FE model (Figure 12 (b) and (c), respectively), it is apparent 
that the lower frequency absorption is enhanced while the higher frequency peak is 
suppressed. This causes the asymmetric peak shown in Figure 12(d), as the two 
absorption peaks are merged by making the dimension of squares close to each other. 
This can partially be explained by the fact that the smaller squares are more active at the 
larger square resonant frequency than vice versa as shown in Figure 13. 



(a) (b) 

Figure 13. FE model solution for 16 and 18 pm squares at the resonant absorption 
frequency of the 18 pm squares (a) and the 16 pm squares (b). The center 
square is 18 pm. Electric (red) and magnetic (blue) fields are shown. Note 
how the 16 pm squares have a larger electric field in (a) than the 18 pm 
squares do in (b) showing that the 16 pm squares assist the 18 pm at their 
resonance, but not vice versa. 


Another configuration for producing two absorption peaks is shown in Eigure 


14(a) where smaller squares are placed around larger ones in a tile-like pattern. This 
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configuration allows for squares with large size differenees, and therefore absorption 
peaks at widely separated frequeneies, to be more densely paeked to increase absorption 
magnitude as demonstrated in Figure 14(b). Measurements showed two distinet 
absorption peaks at frequeneies 4 and 6.5 THz, as depicted in Figure 14(b), with the high 
frequency peak associated with the smaller squares, as expected. Additionally, the higher 
frequency peak has a consistently larger magnitude, most likely due to the higher density 
of 10 pm squares. It is also apparent that the strength of the peak at 6.5 THz drops faster 
compared to the one at 4 THz as the thickness of the dielectric layer is increased. The 
6.5 THz peak seems to be mueh more sensitive to changes in dieleetric thiekness 
compared to single-absorption peak eonfigurations, sueh as in Figure 10(c). 




C) (b) 

Figure 14. Unit eell (a) and absorption charaeteristies (b) of a two peak metamaterial 
layer utilizing a tile eonfiguration to maximize square density. From [23]. 


The number of absorption peaks can be extended from dual-band absorbers by 

introducing additional square sizes. The geometry in Fig. 15 (a) has three different 

squares (11, 14, and 17 pm) and maintains symmetry in two direetions, thus remaining 

polarization-independent at normal incidenee, while also allowing for one of the squares 

(in this case the 14 pm squares) to have a 20% higher density, enhaneing absorption at 

one of the corresponding resonant frequencies. This enhancement can be useful to 

increase absorption selectively if desired. As with the dual-band configuration, three 

distinct peaks associated with each square size can be resolved in the measured 

absorption spectrum in Fig. 15(b). In this case, the middle peak at 5 THz has a higher 
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square density assoeiated with it and an absorption peak of similar magnitude eompared 
to the one eorresponding to the larger squares at 4 THz. The overall absorption was found 
to be lower than that obtained for single-band and dual-band metamaterials, likely due to 
the lower density of partieular square sizes in this eonfiguration. 




Frequency (THz) 

(b) 


Figure 15. Unit cell (a) and measured absorption (b) for a triple band metamaterial 
configuration with square sizes of 11, 14, and 17 pm showing three clear 

peaks. From [23]. 


As metamaterial absorbers have a narrowband absorption peak, Kirchoff s Law 
demands that they also possess a narrowband thermal emission peak [50], [52]. Indeed, 
the emissivity of the metamaterials should follow their absorptivity exactly for thermal 
equilibrium to be upheld. As described earlier, this was tested using the metamaterial 
absorbers as an external source for the FTIR and comparing it to a blackbody. Samples 
A, B, and C consisted of single band absorbers with 17 pm, 13 pm, and 9 pm squares, 
respectively. All samples had 21 pm periodicity and 1.6 pm thick dielectric spacer of 
stoichiometric SiOi. The emission spectrum of a reference emitter (a silicon wafer coated 
with black carbon) was used to extract the emissivity of the metamaterial samples as a 
function of THz frequency. The emissivity spectra, shown in Figure 16 were obtained by 
dividing the sample data by the corresponding data obtained using the blackened emitter 
heated to the same temperature (400 °C). 
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Figure 16. Measured emissivity of the metamaterial samples A, B and C at 400 °C. 

Emissivity exhibits peaks at 4.1, 5.4 and 7.8 THz, respectively. From [24]. 


It is clear that the emission spectra in Figure 16 agree well with the measured 
absorption data in Figure 16 as expected from KirchhofFs law. The emission bandwidth 
at half maximum for each metamaterial structure is approximately 1 THz. 
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The total power emitted into eaeh THz band per unit area (P) can be estimated by 
integrating Planck’s blackbody radiation formula (P(v)) multiplied by the measured 
emissivity (s(v)) within the emission bandwidth as: 


P{v) 


c •’ 


:(v)i 


exp 


hv/ 


'^kT. 


-dv 


-1 


( 3 ) 


where h is Planck’s constant, c is the speed of light in vacuum and k is Boltzmann’s 
constant. The estimated power per unit area, based on Eq. (3), for samples A, B and C 
operated at 400°C is 11, 18 and 36 W/m^, respectively. 



Figure 17. Spectral absorptance of the metamaterial samples A, B and C. Solid lines 
represent the measured absorptance with peaks at 4.1, 5.4 and 7.8 THz, 
respectively. Dashed lines represent the finite element simulation results. 

From [24]. 
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III. METAMATERIALS INTEGRATED INTO THERMAL 

SENSORS 


Note: The first section is extended from a paper published in Optical Materials 
Express [25]. The remaining material in this chapter is contained entirely within a paper 
published in Optics Express [26], and was also presented at the SPIE Photonics West 
conference in February 2013. 

A. THERMAL SENSOR PROPERTIES 

To achieve an increase in performance for thermal sensors, we must first consider 
their performance characteristics. The two most significant parameters are the thermal 
time constant (x) and responsivity (Rth). All thermal sensors have a thermally isolated 
sensitive region which disperses heat to the environment through radiation, convection 
and conduction as illustrated in Figure 18. 


Sensor 

C 

^Leg 


Heat sink 

P 

Figure 18. Diagram showing heat exchange between a thermal sensor and its 

environment. Incoming power P creates a gradient of AT between the sensor 
and the heat sink at temperature T in addition to radiation losses. Convection 
losses are ignored as most sensors operate at low pressures. After [20]. 




Radiation losses are unavoidable and intrinsic to all thermal sensors as KirchofP s 

law [50] demands that any object that absorbs radiation must possess the same emissivity, 

and it is clear that THz metamaterial absorbers are no exception [24]. Convection can be 

controlled through the gas pressure of the surrounding environment [52] and conduction 

through materials and structural design of the sensors. If the thermal losses can be 

approximated as linear with respect to the temperature gradient (AT) between the 
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thermally isolated sensitive area and the surrounding environment, they ean be lumped 
into a single eonstant term ealled the thermal eonductanee (G;/,). Therefore, the heat 
exehange equation for a thermal sensor absorbing radiation of power will be; 


4(0 = G..Ar + Q— (4) 

at 

where Cth is the thermal eapaeitance of the thermally isolated sensitive area. Solving this 
differential equation, we can then arrive at a formula for the responsivity of a detector 
under sinusoidal illumination power [10]; 




AT 


ri 


gJ\ + (oV 


( 5 ) 


where rj is the absorbance, the fraction of incident power {Pg) absorbed by the sensitive 
region of the detector, r is the thermal time constant (r = ) and co is the angular 


modulation frequency of the incident radiation. The graph of Rth vs Gth in Figure 19 
shows that reducing thermal conductance is always advantageous in terms of 
responsivity. 



Figure 19. Responsivity {Rth [K/W]) vs thermal conductivity {Gth [W/K]) in a thermal 
sensor under sinusoidal illumination power. Clearly, regardless of Gth, Rth 
increases with higher absorption (?/), lower frequencies {oS), and smaller Cth 
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Indeed, if thermal eonduetanee were zero the sensor would simply heat up 
indefinitely from ineident radiation. However, real sensors do not have an infinite 
dynamie range and, of eourse, Kirehoff s law ensures that radiation losses eannot be 
avoided [50]. Most thermal sensors are designed sueh that the unavoidable radiation 
losses dominate the thermal eonduetanee to maximize sensitivity [10]. Therefore, 
integrating metamaterials will primarily infiuenee x and rj. Clearly, metamaterials provide 
rj elose to 100 %, but x must be minimized for real time imaging applieations, as thermal 
deteetors typieally operate in a frequeney range where wx « 1 to aehieve the highest 
responsivity. Therefore, if x is large, 0 must be small making real time imaging 
problematie. This means that beeause Gth is already minimized, to minimize x we must 
minimize Cth- The volume of the materials and their thermal eapaeitanee will determine 
the total Cth- The area of the thermal sensor will be determined by the applieation and 
wavelengths of interest; indeed, more sensor area ean result in a higher signal from 
eapturing more ineident power. The thiekness of the thermally sensitive area, however, 
ean be eontrolled through design and should thus be as small as possible. This is the 
primary reason for integrating metamaterials into these sensors, the struetural materials 
are poor THz absorbers by themselves [21], [22], and exeessively thiek struetures will 
perform poorly in a real time imaging setup. 

B. BI-MATERIAL SENSOR CHARACTERISTICS 

As the previous seetion illustrates, the basie physies of integrating metamaterials 
into thermal sensors is similar in that thermal eapaeitanee and absorption are the primary 
eoneerns. Of eourse, eaeh partieular sensor will have its own ehallenges, and the speeifies 
of these ehallenges for bolometers will be explored in the following ehapter. Initially, we 
ehose bi-material sensors to test funetionality metamaterials due to the experienee of my 
eolleagues in their fabrieation and the lab equipment available. 
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Like all thermal deteetors, bi-material sensors respond to a ehange in temperature 
of their absorbing element; in this case, the pixel experiences angular deflection per unit 
power ) and this responsivity can be described as [26]; 

= ( 6 ) 

dP ‘"dT G,,4^^^dT 

Clearly, the expression is simply the temperature response of the sensor multiplied by the 
deflection per unit temperature change (thermomechanical sensitivity). As described 
before, the thermal time constant limits the speed of the sensor. Noise must also be 
examined when evaluating performance; in bi-material sensors noise arises from several 
different sources such as temperature fluctuations, background fluctuations, thermo¬ 
mechanical resonances, illumination source fluctuations and the readout system [53]. The 
first four manifest as fluctuations in the overall sensor deflection, while the readout noise 
depends on the probing mechanism. In a practical sense, the total noise of the complete 
detection system can be described by the noise equivalent power (NEP). For bi-material 
sensors, NEP can be defined as the incident radiant power that produces an angular 
deflection equal to detector’s root mean square (rms) noise [54], [55]. 

Fundamentally, there are two main choices when designing a bi-material sensor; 

materials and configuration. Materials should be fabrication-friendly, exhibit low residual 

stress, have very different thermal expansion coefficients and strong THz absorption. 

Configurations should have a large absorption area, good thermal isolation to increase 

sensitivity, provide a reflective surface for optical readout, and have a small thermal 

capacitance. All of these requirements are intrinsically interdependent making the 

optimization of the final sensor highly dependent on the intended application. 

Nonetheless, the quest to achieve high performance THz bi-material detectors starts with 

dOjdT, defined by the bimetallic effect, and p, which is maximized by using 

metamaterial based absorbers. The sensors here were not the first attempted, in our 

previous work [31] a complex configuration of two fourfold legs metalized in opposite 

sides were used. Although the total length, and thus bi-material effect, was higher than 

that of the current design (see Figure 20(a)), the THz sensitive area was half the size to 

keep the overall pixel sixe around 200 x 200 pm^. The metallization of the front and back 
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sides of alternate legs eaused the strueture to be highly stressed after release. As a result, 
the sensors were bent as much as 30 degrees out of the substrate plane, which is 
problematic for the optical readout. The configuration reported here exhibits several 
advantages, which are highlighted in the next sections. 



Figure 20. Bi-material sensor, (a) 3D view of the THz bi-material sensor with integrated 
metamaterial absorber, fabricated on a Si substrate, (b) Close up of an 
isolated bi-material beam which length is h, and metal and dielectric 
thickness ti and ti, respectively. Azieg and A0 are the linear and angular 
deflection of the beam, respectively. From [26]. 


To increase sensitivity, it is important to optimize the bi-material layer thickness 
to maximize the deflection under increasing temperature. The linear displacement (Azieg) 
of the free tip of a bi-material beam, as depicted in Fig. 19(b), was first quantified by 
Timoshenko [56] and further explored by other groups [53], [57]. If the linear 
displacement is much smaller than the length of the bi-material beam (^ ), which is true 

in our case, the angular deflection due to temperature change {dOjdT) or 
thermomechanical sensitivity can be estimated using; 


dT 


(h 


-2 4 + 6fL + 4fl_ + .‘l . ‘2 -2 




h ^2 


h j 


(a^ - (^2) 


( 7 ) 


where t represents thickness, a is the thermal expansion coefficient and E is the Young’s 
modulus. The indices 1 and 2 are used to represent materials 1 and 2, respectively. Table 
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1 lists some of the most eornmon MEMS materials along with their struetural, thermal 
and eleetrical eharacteristies. 


When the bi-material legs are eonneeted to a freestanding flat absorber (see 
Figure 21(a)), the sensor angular defleetion is approximately equal to A6 (see Figure 
20(b)). The effeet ean be further amplified by adding multifold legs with alternate bi¬ 
material segments [31]. However, such a configuration also magnifies the bending due to 
residual stress after release. 


Material 

Si 

SiNx 

Si02 

A1 

Au 

Young’s Modulus 
E(x 10^ Pa) 
Expansion 

100 

180 

68 

70 

77 

Coefficient 

a(xlO-K-^) 

Thermal 

2.7 

2.1 

0.4 

25 

14.2 

Conductivity 
g (Wm'‘K‘) 

130 

19 

1.4 

237 

296 

Heat Capacity 
c(Jkg‘K-‘) 

750 

691 

703 

900 

129 

Density 

p (x 10- kg m'^) 
Electric 

2330 

2400 

2200 

2700 

19300 

Conductivity 
a (x 10^ S m~‘) 
THz 




10 

37 

refractive index 
* 

n 

3.48-0.11 

2.1-0.0251 

2.0-0.21 

- 

- 


Table 1. Properties of standard MEMS materials[58], [59]. From [26]. 


It is clear from Figure 21 that the Al/Si02 combination produces the highest 
sensitivity with the maximum occurring when the metal thickness is approximately one- 
half of the dielectric thickness. Non-stoichiometric SiNx can provide less stressed layers 
than Si02 [60], however, silicon-rich SiOx can be deposited with much lower stress than 
Si 02 , while preserving most of the thermomechanical and electro-optical properties. 
During sensor fabrication, testing layers of non-stoichiometric SiOx and stochiometric 
Si02 layers with the same thickness were deposited on Si substrates with intrinsic stress 
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on the order of-13 MPa and -140 MPa, respeetively. By seleeting SiOx and A1 (both 
standard MEMS materials), it is possible to maximize dO/dT while simultaneously 
alleviating some of the excessive residual stress related deformation observed in the 
sensors fabricated in [31]. Furthermore, SiOx and A1 exhibit electro-optical properties 
that are suitable for highly efficient metamaterial absorbers, as shown in the previous 
chapter. 


- Al/SiOj -- 

--- Au/SiO^ 

- AI/SiN, -- 

--- Au/SiN_ 

- Al/Si 

--- Au/Si 

• Al/SiOx (FE / Experimental) 




— L 




! Bi-m aterial Leg 






A9 ( 




leg 


(a) 



(b) 


Figure 21. Thermomechanical deflection of the bi-material sensor, (a) Freestanding flat 
THz absorber connected to a bi-material beam, whose length is h, and metal 
and dielectric thickness ti and t 2 ,respectively. Azabs is the total linear 
displacement andzl0 is the angular deflection of the absorber, (b) 
Thermomechanical sensitivity (dO/dT) of the structure of part (a), calculated 
using Equation (4) for all combinations of metal/dielectric of Table 1 where 
ti varies from 10 to 800 nm and t 2 is fixed in 1.1 pm. The circular marker 
shows FE and experimental results for ti = 170 nm. Solid lines are for A1 
combinations while dashed lines are for Au. From [26]. 


C. BI-MATERIAL THZ SENSOR DESIGN 

Bi-material terahertz detectors were designed using a metamaterial structure 
optimized to absorb at 3.8 TFlz radiation form a QCL available to us. Relatively large 
pixel dimensions were chosen to increase the absorption area and simplify the fabrication 
and characterization process. Thermal conductance was intentionally varied among the 
designs while thermal capacitance remained essentially constant (see Table 2). Figure 22 
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shows the geometrical details of three sensors (A, B, and C) with different thermal 
conductances. 
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Figure 22. Structural parameters of the three bi-material THz sensors in pm. (a) Top view 
of sensor A showing all dimensions, (b) Top view of sensors B and C, 
showing the differences in sizes of the thermal insulator anchors, (c) Vertical 
cut of the sensor structure. From [26]. 


The detectors consist of a square metamaterial sensing element in the center 
connected to two symmetrically located rectangular bi-material legs. The entire structure 
is then connected to and thermally isolated from the substrate by folded SiOx anchors 
with varied dimensions as shown in Figure 22. The thickness of the A1 ground plane and 
squares is 100 nm while the bimaterial legs have a 170 nm layer of A1 on the top side. 
The structural SiOx is 1.1 pm thick everywhere. The thermal conductance (G) of all the 
sensors was estimated using the expression; 


Q _ Sth^C 

I 


( 8 ) 


where is the thermal conductivity, is the cross-sectional area and / is the length. 

Since the dimensions of the thermal isolation sections are different, the total thermal 
conductance was estimated by adding the thermal resistance of each section. The 
metalized parts are considered thermal shorts due to their high thermal conductivity 
compared to that of SiOx. The thermal conductance via radiation loss of heat is found to 

be an order of magnitude lower than that via the legs due to low emissivity of A1 and the 
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THz metamaterial that cover most of the sensor surfaces. Heat dissipation due to 
convection is negligible as the sensors typically operate under low pressure (in a vacuum 
chamber). The thermal capacitance was estimated using the expression: 

C = c,,pA^t (9) 

where, is the material thermal capacity, p is the material density. As is the surface 

area and t is the layer thickness. The thermal capacitance of the sensor is the sum of 
thermal capacitances of the SiOx and A1 layers. The material parameters used for the 
calculations are given in Table 1. The time constant (x = C/G) was also estimated for each 
sensor configuration and listed in Table 2 in addition to other parameters. 


Sensor 


Sensor A 



SensorB 



Sensor C 


Property 

Anal 

FE 

Exp. 

Anal 

FE 

Exp. 

Anal 

FE 

Exp. 

Absorptance 

n 

- 

0.96 

0.95 

- 

0.96 

0.95 

- 

0.96 

0.95 

Thermal 

Conductance 

1.6 

1.7 


2.2 

2.1 


9.3 

8.5 


G (x 10- WK^) 
Thermal 
Capacitance 

11.1 

12 


10.7 

12.5 


9.8 

11.9 


C (x 10- JK') 
Time Constant 

0.68 

0.7 

0.8 

0.47 

0.6 

0.5 

0.1 

0.14 

0.3 

t(s) 

Thermomechanic 
al Sensitivity 

0.15 

0.19 

0.2 

0.15 

0.2 

0.2 

0.15 

0.25 

0.2 

dO/dT (deg K^) 
Responsivity 
dO/dP (x 10'’ deg 

0.95 

1.1 

1.2 

0.65 

0.9 

0.8 

0.15 

0.25 

0.2 

W‘) 

Noise Equivalent 
Power 

(due to incident 

0.00 


8.6 

0.00 

f. 


13 

0.01 

A 


45 

power) 

NEP (x 10- W) 




o 



H 




Table 2. THz bi-material sensor analytical, numerical, and experimental parameters. 

From [26]. 
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The deformation of the sensor with inereasing temperature was analyzed using the 
COMSOL heat transfer module, which allows a uniformly distributed heat flux boundary 
to be placed at the absorber to emulate the incoming THz power. The anchor attachments 
to the substrate are fixed and set at a constant temperature to represent the heat sink. All 
other boundaries are thermally insulated from the surroundings and free to move. The 
program computes the heat transfer equation at each mesh point allowing the retrieval of 
several parameters, such as temperature distribution, thermal deformation, etc. For steady 
state simulations the total incoming heat flux was conveniently set as 1 pW, therefore the 
thermal deformation and temperature distribution can be directly read “per unit pW.” The 
angular deformation can be directly obtained by the displacement of the free edges of 
absorber and hence d0/dT can be estimated using the temperature difference between the 
absorber and heat sink. Also, the responsivity (dO/dP) of the sensors can be obtained 
using the maximum deformation (steady state) and the incident heat flux (1 pW). 
Furthermore, thermal conductance can be estimated using Equation 8. 

Time domain simulations were performed to obtain the transient response of the 
sensor to a pulsed heat flux allowing the retrieval of the time constant of the sensors. 
Using the obtained time constant and thermal conductance, we estimated the thermal 
capacitance of the sensors. The calculated and simulated parameters, using the material 
properties of Table 1, are listed in Table 2 and, in general, show good agreement. Notice 
that the thermal capacitance values obtained by FE simulations show a small discrepancy 
as they increase with decreasing sensor mass. This is most likely due to the time constant 
estimation, which is more susceptible to errors as it decreases. Eigure 24(a) and 24(b) 
show the deformation plots obtained by EE simulation, where the z-axes are scaled up 20 
times for visual purposes. The surface color scale indicates the temperature distribution 
and it is the same for all sensors. It can be seen in Eigure 24 that sensor A deflects more 
compared to sensors B and C under the same incident power (1 pW) primarily due to 
lower thermal conductance. Eigure 24(d) shows the time domain simulations where the 
sensors are submitted to a step excitation (black solid line) of 1 pW for duration of 8 
seconds. The vertical axes show temperature on the left side and angular deflection on the 
right side. 
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Figure 23. FE simulations showing the deformation plots of sensor A (a), B (b) and C (c) 
under a eonstant 1 pW heat flux. The z-axes are sealed up 20 times for visual 
purposes. The surface colors indicate the temperature distribution according 
to the color bar on the left, (d) Time domain simulation of all three sensors 
under a 1 pW step excitation (black line). Temperature change and angular 
displacement are shown on the left and on the right, respectively. From [26]. 


Noise sources intrinsic to the detectors were also considered and an analysis 
similar to that in [53] was performed to determine the NEP. The expressions given by 
Equations 10 and 11 were adapted from [53] to reflect angular deflection fluctuations. 
The primary noise sources in thermal detectors are temperature fluctuation, background 
fluctuation and thermomechanical noises. The spontaneous fluctuation in angular 
deflection (deg) of the absorbers caused by temperature fluctuations is given by; 

( 10 ) 

\ / ^ 

where T is the sensor temperature, ke is the Boltzmann constant, G is the total thermal 
conductance and B is the bandwidth, which can be set to unity. The background 
fluctuation noise can be obtained by replacing the total thermal conductance in Equation 
10 by thermal conductance via radiation loss of heat. Flowever, this is much smaller than 
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the thermal conduetance via the legs and its contribution to noise can be neglected. The 
angular deflection (deg) due to thermomechanical noise, knowing that the detector 
operating frequency is much slower than the mechanical resonances (few kHz), is given 
by: 



360 Ukj^GB 

Tzl^ ]/ Qko), 


( 11 ) 


where Q is the quality factor, k is the stiffness and 0o is the resonant angular frequency of 
the mechanical structure. Using the eigenfrequecy solver in the COMSOL structural 
mechanics module, the first resonant frequency and stiffness of all the sensors were 
estimated and found to have values 3.5, 4.0 and 6.0 kHz and 0.02, 0.025 and 0.04 Nm~' 
for sensors A, B and C, respectively. Typical Q values for similar structures lie between 
100 and 1000 in vacuum [53]. The noise was estimated and as expected, the dominant 
source is the temperature fluctuation in the detector. The total noise intrinsic to the 
sensors was estimated to be 5.0, 4.0 and 2.0 pdeg. The NEP values of the three sensors 
were calculated by dividing the fluctuations due to the noise by their respective 
responsivities, and are listed in Table 2. 


D. FABRICATION AND CHARACTERIZATION 

The sensors were fabricated using standard micromachining technology. First, a 
100 nm thick aluminum (Al) film was deposited on a 300 pm thick silicon (Si) substrate 
by e-beam evaporation. Then, the Al layer was patterned and wet etched to form the 
absorber ground plane. Next, a 1.1 pm thick SiOx layer was deposited using plasma 
enhanced chemical vapor deposition (PECVD) at 300 °C, followed by another 100 nm 
thick Al film. The second Al layer was then patterned and plasma etched to define the 
absorber metamaterial squares. Then a 170 nm thick Al layer was deposited, patterned 
and lifted off to form the bi-material legs. The sensor structure was then created by 
reactive ion etching of the SiOx layer. Finally, the structures were released through 
backside trenching using the Bosch etch process. Circular openings were chosen to 
ensure release of the structure and to help refine the Bosch etch recipe. Figure 24(a) and 
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24(c) show the optical profile of the fabricated sensor A along with micrographs of the 
three sensors and an array of sensors of type A. 



Figure 24. Fabricated THz bi-material sensors, (a) 3D optical profile of sensor A (the 
aspect ratio is preserved), (b) Micrograph of an array of sensor A. (c) 2D 
profile taken along the bimaterial legs direction (y-profile) with the 
processing direction (z-profde) scale exaggerated to show the residual 
deflection of the legs (red line) and absorber (blue line), (c) Micrographs 
showing the top view of sensors A, B and C. From [26]. 


The 3D profile in Figure 24(a) is an actual view of sensor A while the 2D profile 
in Figure 24(c), taken along the bi-material legs direction (y-profile), is exaggerated in 
the processing direction (z-profile) to show the residual deflection of the legs and 
absorber of the sensor. The measured residual deflection of the absorber is approximately 
6° for the sensors A and B and 8° for sensor C. It is easy to observe that the absorber is 
almost flat due compensation of stresses from the aluminum layers in both sides of the 


43 
















SiOx layer. Due to the defleetion of the sensors, mierographs shown in Figure 24(d) are 
not eompletely focused across the surface. In addition to the sensors, the fabricated wafer 
contains an area of 10 x 10 mm filled with the same metamaterial structure used in the 
sensors. This is to allow accurate measurement of the absorption characteristics of 
metamaterial used in the sensors. The absorptance of the metamaterial film was measured 
as described in chapter II and compared with the QCL emission characteristics as shown 
in Figure 25(a). A good match between the absorptance peak position of the metamaterial 
and the 3.8 THz QCL emission frequency was achieved with nearly 90% absorptance. 
This assured that the sensors absorbed our QCL emission with high efficiency. 



Figure 25. (a) Measurement of the absorptance spectra of the THz sensors metamaterial 
structure (blue line) compared with the QCL normalized emission (read line), 
(b) Measured angular deflection (markers) upon temperature chance. Notice 
that the effect is linear and almost indistinguishable among the sensors, 
resulting in approximately 0.2 deg/K thermomechanical sensitivity. From 

[26]. 


Next, the thermal response of the sensor {dO/dT) was measured. The temperature 
gradient in the bi-material section of the leg was estimated to be less than 5% of that 
between the central absorbing element and the substrate. Thus, the bi-material section of 
the leg can be treated as thermally shorted allowing the measurement of the thermal 
response by uniformly heating the sensor. The measurement was performed by attaching 
the sensor to a fiat resistive heating element and sweeping the temperature from 303 to 
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313 K. The reflection of a laser diode beam from the backside of the sensor’s ground 
plane was projected on a screen and the angular deflection of the sensor was determined. 
Angular deflections from the three sensors are shown in Figure 25(b) with different 
markers. The deflections are almost indistinguishable because the detectors have the 
same bi-material leg dimensions. The solid line is a linear fit, showing that the response 
in this temperature range is linear and approximately 0.2 deg/K, which is slightly higher 
than the estimated values (see Table 2 and Figure 21). Figure 21 shows that the thermal 
response of the sensors can be further increased by, for example, decreasing the SiOx 
thickness or increasing the A1 thickness on the legs. Test structures fabricated in parallel 
with these sensors showed that increasing the A1 thickness on the legs also increases the 
residual stress. Decreasing the dielectric thickness has a similar effect in addition to 
reducing absorber efficiency [30]. Further tuning of the fabrication process is necessary 
to reduce the residual stress on the bi-material legs to decrease the initial bending as 
depicted in Figure 7(a). 

Subsequently, the sensors were placed in a vacuum chamber and operated at a 
pressure of approximately 0.03 mTorr to minimize the heat loss by convection [52]. The 
QCL was kept inside a cryostat and operated at around 15 K. The divergent THz beam 
passed through the cryostat Tsurupica window and the radiation was focused by a 40 mm 
polyethylene lens onto the sensors. Both Tsurupica and polyethylene exhibit reasonable 
transmission (-65%) in the THz range. The QCL was operated in pulsed mode with the 
pulse width fixed at 5 ps and a variable pulse rate to control the output power. The 
deflection of the sensor was measured using the same procedure described earlier for a 
set of QCL pulse rates ranging from nearly zero to 5 kHz. The absolute power that 
reaches the sensors (incident power) is estimated using the responsivity {dO/dP) in 
Equation 6 along with the calculated thermal conductance and measured absorptance. 
Note that the QCL switching frequency and duty factor must be taken into account since 
the sensors can only respond to the average power. Figure 26(a) shows the measured 
angular deflection versus incident power for all three sensors (colored markers). For all of 
the sensors, the responsivity values estimated analytically tend to be lower than that of 
the FE and experimental values mainly due to the underestimation of bimetallic effect 
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(dO/dT), by Equation 7 (see Figure 21). As expected from Equation 6, responsivity of the 
sensors was found to decrease with increasing thermal conductance. 



Figure 26. Responsivity and NEP measurements, (a) Measured angular deflection per 
varying incident power for all three sensors (colored markers) Notice that 
responsivity increases as thermal conductance decreases, (b) Measured output 
voltage of the PSD for sensor A by gating the QCL output at 200 mHz. The 
power incident in the detector is shown on the right vertical axis. From [26]. 


To determine NEP, a position-sensing detector (PSD) was added to the 
experimental setup to read the deflection at low power levels. The NEP was then 
measured for each detector and listed in Table 2. Figure 26(b) shows measured output 
voltage of the PSD for sensor A by gating the QCE output at 200 mHz. It is important to 
highlight that the measurements include the effects of QCE power fluctuations and 
optical readout noise, not considered in the theoretical estimations discussed earlier. The 
difference between the measured values (3 orders of magnitude higher) and the estimated 
ones (Table 2) can be attributed primarily to the readout noise. The QCE power 
fluctuations [61] do not seem to contribute to the observed noise since the noise floor 
when the QCL is off, shown in Figure 26(b), is similar to the noise observed when the 
QCL is on. As expected, NEP increases from sensor A to C due to decrease in 
responsivity. The measured NEP values, including the readout noise and the intrinsic 
noise of the sensor, can be translated into minimum detectable temperature difference on 
the sensor, found to be approximately 50 mK for all three sensors. This value is similar to 
those of bi-material sensors operating in the IR range [53], [58]. 
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The time domain response was also measured using the PSD and the results for 
the three sensors are shown in Figure 27(a) under the same ineident power with the QCL 
gated at 500 mHz. As observed in Figure 26, sensor A is more sensitive, whieh agrees 
with the predictions and previous measurements. Since the sensors have the same p, the 
same absorbing area, same materials, the same dO/dT, and nearly the same thermal 
capacitance, speed and responsivity are completely controlled by the thermal 
conductance, which depends on the anchor geometry. 




Figure 27. Time and frequency domain measurements, (a) Time responses of sensors A, 
B and C measured under the same incident power with the QCL gated at 1 
Hz. Noticed that sensor A is more sensitive and slower, which agrees with the 
predictions and previous measurements, (b) Normalized frequency responses 
for the three sensors (colored lines). The time constants were retrieved by 
taking the inverse of the 3 dB frequencies that are 1.2, 2.1 and 3.2 rad/s for 
sensor A, B and C, respectively. From [26]. 


The time constant of the sensors was determined by sweeping the QCL gating 
frequency from 50 mHz to 30 Hz and recording the PSD peak to peak voltage. The 
normalized frequency responses for the three sensors are shown in Figure 10(b). The time 
constants were retrieved from Figure 27(b) by taking the inverse of the 3 dB frequency 
and included in Table 2. In general, the measured time constants agree well with the FE 
estimations, while the analytical approach underestimates this parameter. The sensors are, 
apparently, slow for real time imaging; however, pixel size can be reduced to increase the 

speed of operation. A tradeoff between speed and sensitivity has to be made according to 
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specific applications. Although the fabricated sensor arrays do not have high spatial 
resolution, their imaging capabilities were probed by a CCD camera with coaxial 
illumination as schematically illustrated in Figure 28(a) [62], [63]. An image of the QCL 
beam, gated at 500 mHz, is shown in Figure 28(c) compared with the same image 
obtained using a commercial IR microbolometer camera with THz optics [64] (Figure 
27(b)). The focal plane array of the IR camera has 30 pm pitch of and can resolve the 
rings associated with the QCL beam. Our sensor array, on the other hand, has a 430 pm 
pitch and cannot resolve the rings; nevertheless, it gives a raw image that clearly shows 
where the energy is concentrated and the circular shape of the THz beam. Finally, it is 
possible to optimize THz bi-material sensor configuration, size, fabrication processes and 
readout to achieve real time imaging. 
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(b) (c) 


Figure 28. QCL beam imaging, (a) Optical readout used to record videos [46] and the 
snap shot showed in part (c). The images were recorded using background 
subtraction to suppress to the effects of the residual stress of the sensors, (b) 
Image obtained using a 30 pm pitch commercial IR microbolometer camera 
with THz optics, (c) Image (Media I) of the same QCL beam, gated at 500 
mHz, obtained using an array of sensor A with 430 pm pitch using the 
readout depicted in part (a). Notice that since the pitch of our sensor are one 
order of magnitude higher than the IR camera, it cannot resolve the rings 
associated with the QCL beam, showed in part (b). From [26]. 
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IV. INTEGRATING METAMATERIALS INTO 
MICROBOLOMETER PIXELS 


Note: Other than the last section, this chapter is published practically unaltered 
in Optical Materials Express [25]. 

A. CHALLENGES TO BOLOMETER INTEGRATION 

As the previous chapters have described, metamaterials can make effective thin 
THz absorbers for thermal sensors. Bi-material sensors demonstrated in Chapter III the 
metamaterial THz absorber and bi-material legs are physically separated and do not 
interfere with each other’s function. In the case of microbolometers, a temperature 
sensitive layer must be included in the central absorbing area, or microbridge which can 
affect the efficiency of the metamaterial absorber. This resistive bolometric layer can be 
made out of several different materials, such as titanium (Ti) and vanadium oxides 
(VaOx). In Chapter III, it was shown that the absorbing layers should be as thin as 
possible to reduce thermal capacitance {Cth), and thus the thinnest choice would be one in 
which the bolometric layer is simply integrated into the middle of the dielectric spacer as 
in Figure 29 or serve as the ground plane of the metamaterial. Unfortunately, the latter 
option leads to a significant reduction in absorption as THz will easily penetrate high 
resistivity bolometric materials and thus defeats the purpose. Therefore, the former option 
was investigated, although some other possible configurations will be examined. 

B. FABRICATION AND MEASUREMENT 

The approach used in this work involves a metamaterial structure with a 
bolometric layer embedded in the dielectric layer as schematically illustrated in Figure 
29. The aim of the work in this chapter is to determine the optimal location of the 
bolometric layer within the dielectric to achieve the highest possible THz absorption. 
Three wafers containing metamaterial absorbers with Ti as the bolometric layer located in 
the dielectric spacer at varying depths along with one without a bolometric layer were 
fabricated using standard MEMS microfabrication techniques, similar to those used in 
Chapter II. The fabrication sequence includes deposition of an 85 nm thick layer of A1 on 
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4” Si substrates using e-beam evaporation. Following this, a layer of low stress 
nonstoehiometric silicon-rieh silieon oxide (SiOx) with different thieknesses was 
deposited using plasma enhanced chemical vapor deposition (PECVD). 




Figure 29. Schematic (a) and top view micrograph (b) of a metamaterial absorber with a 
bolometric layer embedded in the dielectric. The square size and pitch were 
selected to have an absorption peak within the 1-10 THz band. From [25]. 

Wafers A-D possessed 1270 nm, 765 nm, 297 nm, and 1590 nm of SiOx, 
respectively, confirmed through optical interferometry. Then, a Ti layer was sputter 
deposited using Ar plasma onto wafers A-C and a Si test wafer with a SiNx isolation layer 
to measure conductivity. The measured thickness of Ti was found to be about 10 nm and 
confirmed to have a sheet resistance of about 220 f2/sq, which corresponds to 
conductivity of 4.45 X 10^ S/m. Following the deposition of Ti, an additional SiOx layer 
with thicknesses of 297 nm, 765 nm, and 1270 nm were deposited using PECVD on 
wafers A-C, respectively. This resulted in nearly equal total spacer thicknesses, including 
Ti layer, of 1577 nm, 1540 nm, 1577 nm, and 1590 nm for wafers A-D. A second 85 nm 
thick layer of A1 was then deposited onto wafers A-D using e-beam evaporation and was 
then patterned using photolithography and sputter-etched with Ar plasma to form an 
arrays of squares as in Chapter II, with the absorption characteristics of the 15 pm side 
and 21 pm pitch configuration shown here. 
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As described in Chapter II, square size mainly influences the resonant absorption 
frequency, while dielectric thickness and periodicity strongly affects absorption 
magnitude. This particular configuration was chosen for its resonant frequency in the 1- 
10 THz range, and other configurations with varying resonant frequencies can be found 
in Chapter II. All wafers were intentionally over-etched approximately 150 nm into the 
SiOx to assure complete insulation between square elements. The composition of layers 
of wafers A-D are summarized in Table 3. To determine the optical properties of these 
metamaterial absorbing layers in the THz spectral range, reflection measurements were 
taken using a Thermo-Nicolet Nexus 870 Fourier Transform Infrared Spectrometer 
(FTIR) with a globar source and PIKE Technologies MappIR accessory. An Au-coated Si 
wafer was used to establish the background for the reflectance measurements. The two 85 
nm thick A1 layers block any significant transmission due to the 40 nm skin depth at 
around 4.7 THz, resulting in over two orders of magnitude reduction of power due to 
exponential attenuation. The periodicities used are less than one third of the free space 
wavelength, resulting in negligible higher order scattering [43]. As before, this allows 
A+T + R =1 to be simplified to A ~1- R . The resulting absorption spectra for samples 
A-C are plotted in Figure 30(a). 


Wafer 

Top SiOx 

Bottom SiOx 

Total 

Spacer^ 

Max. 

Absorption^ 

A 

297 

1270 

1577 

62% 

B 

765 

765 

1540 

88% 

C 

1270 

297 

1577 

94% 

D 

— 

— 

1590 

98% 

“Includes 10 nm of Ti. 

*For 15 pm squares with 21 pm periodicity 


Table 3. Thickness of the dielectric spacer layers (in nm) of each sample. From [25]. 
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(a) 


(b) 


Figure 30. (a) Experimental (solid) and FE model (dashed) absorption speetra for 

samples A (blaek), B (red), and C (green), (b) Theoretieal absorption curves 
for sample A for a set of sheet resistivites in the bolometric layer. Square size 
is 15 pm with 21 pm periodicity for both plots. From [25]. 


C. RESULTS AND DISCUSSION 

It is clear that in Figure 30(a), as the Ti conductive layer gets closer to the 
patterned top A1 layer, the THz absorption of the metamaterials is significantly reduced 
and slightly red shifted. This is expected conceptually as the Ti layer would function as a 
poor ground plane if the lower SiOx and A1 layers were removed. The redshift in the 
resonant absorption frequency is larger than the frequency shifts observed in Chapter II 
for thinner dielectric thicknesses, suggesting that a simple reduction of dielectric 
thicknesses will not explain these absorption spectra. To model the THz absorption 
characteristics of the metamaterial film with an embedded bolometric layer, finite 
element (FE) simulations were performed using the RF module of COMSOF 
multiphysics software explained in greater detail in Chapter II. The periodic geometry of 
the metamaterial absorbing film allows for a single unit cell to be modeled with periodic 
boundary conditions. Internal ports or scattering boundary conditions allow incoming 
plane waves at specific frequencies to enter the model. The THz absorption can be 
calculated by subtracting off reflection and transmission coefficients calculated from the 
S parameters. Alternatively, resistive heating can be used to calculate absorption. This 
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method is especially useful for calculating the absorption of elements that are not 
thermally linked, such as a ground plane disconnected from a bolometer microbridge. 
The aluminum layers were modeled using a conductivity of 1 x 10^ S/m, confirmed 
through four-point probe and stylus profilometer measurements. SiOx was modeled as a 
dielectric with a complex index of refraction of 1.95 + 0.025i, adjusted from the SiOi 
values used in Chapter II based on the experimental results therein. Ti was modeled as 
having a sheet resistance of 220 H/sq based on four-point probe measurements. The exact 
thickness of the Ti is irrelevant as long as the sheet resistance is well known [48] and was 
estimated at 10 nm from previous depositions. As in Chapter II, the over-etching of SiOx 
is included in the models, producing a slight blue shift in the resonant absorption 
frequency due to an effectively smaller index of refraction. Figure 30(a) shows the 
comparison of simulated and measured data for the three samples. The simulated data is 
able to capture the effect of the location of Ti layer on the absorption reasonably well. E- 
field magnitude plots for these simulations are shown in Figure 31, all using the same 
color scale for E-field intensity. 



ik^e 


(b) 


(c) 


Figure 31. E-field magnitude plots from COMSOL EE simulations of 15 pm squares 
with 21 pm pitch at 4.7 THz for samples A (a), B (b), and C (c). All plots use 
the same color scale, with warmer colors indicating more intense E-fields. 
Material types of different regions are labeled in (a). E-field, H-field, and 
propagation vectors of the incident THz waves are shown on (b). From [25]. 
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It is clear that the E-fields do penetrate the Ti layer but not the ground plane in all 
eases, and that the Ti is indeed aeting as a poor ground plane. The E-fields also beeome 
less intense as the Ti layer nears the square array, indieating that it is negatively affeeting 
the resonanee. In addition, absorption of sample A is also simulated for a set of 
eonduetivities of the bolometrie layer as shown in Eigure 30(b). It ean be seen that, if the 
eonduetivity of this layer is redueed, the influenee of the layer on the absorption of the 
metamaterial deereases, as expeeted. Sinee typieal high TCR bolometrie layers have 
sheet resistanees in the range, we should expeet that a spaeer of 300 nm of SiOx 
between the bolometrie layer and the patterned A1 layer will be suffieient for aehieving 
high THz absorption. The results indieate that the two most important parameters that 
affeet the absorption peak height are the eonduetivity and position of the bolometrie layer 
between the patterned A1 layer and ground plane of the metamaterial. The data in Eigure 
30 also seems to show that within the range of parameters used, there is some 
equivalenee in the effeets of sheet resistanee and position of a thin bolometrie layer on 
the absorption speetrum. Qualitatively, less eonduetive layers should allow for greater 
penetration of E-fields and thus more of the desirable resonanee between the square array 
and the ground plane. Also, a high eonduetivity layer near the ground plane or a low 
eonduetivity layer near the square array should both eonverge to the original absorption 
speetrum. 

D. OTHER POSSIBLE CONFIGURATIONS 

Now that we have established that the impaet on absorption is redueed as the 
sheet eonduetivity is redueed and as it is moved eloser to the ground plane, we ean 
evaluate other possible eonfigurations investigated by other researehers as depieted in 
Eigure 32. Maier et al. fabrieated a strueture analogous to Eigure 32(a) [65] and Eigure 
32(e) [32] for the mid-IR range (3-8 pm) range. They note that the bolometrie layer, in 
their ease polysilieon (poly-Si), interfered with their metamaterial signifieantly for the 
eonfiguration in Eigure 32(e). Certainly, the eonfiguration in Eigure 32(a) eliminates 
interferenee between metamaterial and the bolometrie layer sinee they are optieally 
separated sinee the ground plane permitting no signifieant THz transmission. However, 

when eomparing THz metamaterials to mid-IR metamaterials, an important distinetion 
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emerges. Note that in Chapter II, the optimal dielectric thickness for high THz absorption 
was approximately 1 pm, while in the case of Maier et al. the optimal dielectric thickness 
can be as thin as 250 nm [65] also seen in other similar work [42], Thus, the 
configuration in Figure 32(a) is not suitable for THz applications due to increase of 
thermal capacitance when extra layers added for separating the bolometric layer. 



(a) (b) (c) 


Figure 32. Three possible configurations for integrating a conductive layer into a THz 
absorbing metamaterial using a metal, bolometric layer, and a dielectric with 
the bolometric layer below (a) and within (c) the metamaterial layer. 
Configuration (c) thermally disconnects the ground plane. The configuration 
we used is shown in (b) with Ti. 

Of course, thickness can be reduced with some loss in absorption if required and 
controlled over-etching can eliminate some of the dielectric spacer. However, these 
methods will still reduce the performance of the microbolometer pixel; this is why the 
configuration in Figure 32(b) was emphasized for THz metamaterials in addition to the 
fact that, as with the bi-material pixels, the metal layers on both sides of the dielectric 
will compensate for each other in terms of residual stress. The structure depicted in 
Figure 32(c) does present an interesting opportunity, albeit with some fabrication 
challenges. In this configuration, the thickness of the metamaterial absorbing layer 
includes the air gap thus allowing the thermal capacitance of the microbridge to be 
relatively small. The issue here is shown in Figure 5 in Chapter II, the FE model indicates 
that the ground plane has a non-negligible power dissipation of over 20% of the absorbed 
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THz at the resonant frequency. As the dissipation in any material will depend on the 
electric field and conductivity or complex refractive index of that material (assuming 
magnetic losses are negligible), it should be possible to effectively control where the 
absorption occurs. Indeed, if the ground plane were constructed out of a material with a 
high THz conductivity, dissipation in that layer would be less than the other layers, 
minimizing this problem. Therefore, while the structure in Figure 32(b) is likely the 
safest configuration to fabricate microbolometers with high performance. Figure 32(c) 
would allow for an even smaller thermal capacitance. 
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V. CONCLUSION 


The research presented here demonstrates the ability of metamaterials to enhance 
thermal sensors operating in the THz region of the spectrum by improving their 
absorption characteristics. Finite element models provided key insight leading to 
fabricated metamaterial absorbers and bi-material sensors that performed well in the 3-8 
THz range. 

The THz metamaterial absorbers, composed of a metal layer patterned into a 
periodic array of squares separated from a ground plane of continuous metal by a 
dielectric spacer were found to have nearly 100% absorption at a specific THz frequency. 
The peak absorption frequency was inversely proportional to the dimension of the 
squares. In addition, metamaterial THz absorbing layers with one, two, and three 
absorption peaks we fabricated using complex unit cell structures instead of single square 
unit cells. The use of standard microfabrication materials, A1 and SiOx for the fabrication 
of metamaterials, allows them to be incorporated into the construction of thermal sensors 
such as microbolometers or bi-material sensors. The results indicate that these 
metamaterial absorbers can be used as absorbing layers in thermal imagers paired with 
narrow-band sources to significantly enhance their sensitivity through resonant 
absorption at the emitter frequency. In addition, dual- and triple-band configurations can 
be used with multiple lasers or broader sources, as desired. 

The metamaterials were successfully integrated into bi-material pixels. Sensor 
materials and configurations were chosen in order to maximize responsivity. The 
combination of favorable thermal, mechanical and optical properties of the MEMS 
fabrication-friendly materials SiOx and A1 were advantageous. Analytical and FE models 
were used to predict the performance of the sensors. A highly efficient metamaterial 
structure was developed to provide near 100% absorption at 3.8 THz while 
simultaneously serving as a structural layer and providing access for external optical 
readout. The fabricated sensors showed responsivity values as high as 1.2 deg/pW and 
time constants as low as 200 ms, depending on the configuration. Minimum detectable 

power on the order of 10 nW was observed, demonstrating that the sensors can operate 

59 



with low-power THz sources. Although the sensors were not optimized for imaging, the 
use of an external optical readout allowed us to obtain raw images of the QCL beam 
indicating the potential of these detectors to be further optimized for use in focal plane 
arrays for real time THz imaging. The successful design, fabrication, and characterization 
of bi-material MEMS sensors, using metamaterial absorbers demonstrates the potential of 
this approach for THz imaging. 

The final phase of the research, fabrication of THz-absorbing Al/SiOx/Al 
metamaterial layers with a thin Ti conductive layer inserted into the dielectric spacer to 
be used as bolometric layer was explored. The introduction of Ti reduced the magnitude 
of the absorption peak when placed near the patterned A1 layer in addition to a small 
redshift. The measured absorption spectra are in reasonably good agreement with those of 
the simulations. The results indicate that the bolometric layer should be placed near the 
ground plane for achieving maximum absorption. FE simulations of the effects of the 
bolometric layer on absorption suggest weak position sensitivity when the sheet 
resistance of the bolometric layer exceeds about 1 kH/sq. The findings indicate that 
metamaterial THz absorbers should be able to improve the performance of 
microbolometer cameras in the THz range without compromising their thermal time 
constants. 

A. FUTURE WORK 

The knowledge gained in this work can prove useful for the fabrication of 
microbolometer pixels with integrated THz metamaterial absorbers. Therefore, the next 
step in this research should be the fabrication and characterization of microbolometer 
pixels using the optimized metamaterials with embedded bolometric layers for THz 
imaging. To detect a change in resistance in a microbolometer pixel, a Wheatstone bridge 
is employed in readout systems [10] as illustrated in Figure 33. To ensure the load 
resistors are of similar resistance, they can simply be other microbolometer pixels that are 
at ambient temperature. Self-heating is an inevitable issue with a microbolometer 
readout, given that any power dissipated while probing the Wheatstone bridge will 
change the resistance of the microbolometer pixels. Therefore, a possible configuration 
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illustrated in Figure 33 has two thermally shorted pixels, one pixel that is blind to THz by 
joining the squares together (but maintains the same thermal eapaeitanee as the sensitive 
pixel). 



Figure 33. Wheatstone bridge configuration for THz metamaterial microbolometer 
characterization. The thermally shorted pixels are not released from the 
substrate and the blind pixel has a single metal square instead of the patterned 
layer and thus will not absorb THz efficiently, but will have essentially the 
same thermal properties as the sensitive pixel. Several Wheatstone bridges 
are fabricated together, their common ground is not visible. 


This configuration will automatically cancel much of the self-heating, allowing 
for an easier readout. Additionally, a separate configuration with three thermally shorted 
pixels and one sensitive pixel can be used to evaluate the self-heating during readout. In 
commercial systems, an array can be read out pixel by pixel with the same Wheatstone 
bridge, but such readout systems require more fabrication steps to incorporate the 
necessary electronics. 
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